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ABSTRACT 

We report on a survey of the inner part of the Galactic Plane in very high en- 
ergy 7-rays with the H.E.S.S. Cherenkov telescope system. The Galactic Plane 
between ±30° in longitude and ±3° in latitude relative to the Galactic Gentre 
was observed in 500 pointings for a total of 230 hours, reaching an average flux 
sensitivity of 2% of the Crab Nebula at energies above 200 GeV. Fourteen pre- 
viously unknown sources were detected at a significance level greater than 4(T 
after accounting for all trials involved in the search. Initial results on the eight 
most significant of these sources were already reported elsewhere (Aharonian et 
al. 2005a). Here we present detailed spectral and morphological information for 
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all the new sources, along with a discussion on possible counterparts in other 
wavelength bands. The distribution in Galactic latitude of the detected sources 
appears to be consistent with a scale height in the Galactic disk for the par- 
ent population smaller than 100 pc, consistent with expectations for supernova 
remnants and/or pulsar wind nebulae. 

Subject headings: catalogs - surveys - gamma-rays: observations - ISM: super- 
nova remnants - stars: pulsars: general - stars: supernovae: general 

1. Introduction 

Very high energy (VHE, defined here as E > 10^^ eV) 7-rays provide the most di- 
rect view currently available into extreme environments in the local universe. As 7-rays 
are produced by interactions of relativistic particles they represent a good probe of non- 
thermal astrophysical processes. Relativistic electrons can produce 7-rays by non-thermal 
bremsstrahlung or by inverse Compton scattering on background radiation fields, whereas 
protons (and atomic nuclei) can generate 7-rays via the decay of 7r°s produced in hadronic 
interactions with ambient material. VHE 7-radiation is therefore expected from the accel- 
eration sites of particles with energies above 1 TeV (lO"*^^ eV). In particular, 7-rays can be 
seen as a probe of sources of nuclconic cosmic rays in our Galaxy and may thus provide part 
of the solution to the century old puzzle of the origin of this radiation. 

Observations of synchrotron emission in the X-ray band provide strong indirect evidence 
for the existence of ~ 100 TeV electrons in shell-type supernova remnants (SNR), see for 
example Koyama et al. (1995), and in pulsar wind nebulae (PWN), see e.g. Weisskopf et 
al. (2000). At TeV energies, shell-type SNRs and pulsars are established to be sources of 
VHE 7-rays as shown by highly significant detections of the Crab Nebula by the Whipple 
collaboration (Weekes et al. 1989) or the detection of the Supernova remnant RX J1713.7- 
3946 by the CANGAROO (Enomoto et al. 2002) and the H.E.S.S. collaborations (Aharonian 
et al. 2004a). However, despite the fact that SNRs are the best candidates for the acceleration 
of hadronic cosmic rays up to the so-called knee in the cosmic ray particle spectrum at 
~ 10^^ eV, there is as yet no generally accepted evidence for the existence of energetic nuclei 
(as opposed to electrons) in these objects. SNRs are favoured as the site of cosmic-ray 
acceleration for two principal reasons; first, that a theoretically well established acceleration 
mechanism exists in the form of diffusive shock acceleration at SNR shock fronts, for reviews 
see e.g. Blandford & Eichler (1987); Malkov & Drury (2001), and secondly, that the Galactic 
supernovae and the resulting SNRs are the only known potential sources which can provide 
the necessary amount of energy (Ginzburg & Syrovatskii 1964). 
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To explain the observed energy density of cosmic rays, supernova explosions must con- 
vert kinetic energy to energy of relativistic particles with a conversion efficiency of roughly 
0.1 — 0.2. Theoretical expectations for the 7-ray -emission of SNRs arising from hadronic 
interactions are given by Drury, Aharonian & Volk (1994) as 



where d denotes the distance to the object, n the density of the circumstellar medium 
assumed to be uniform, and E'sn the total kinetic energy released by the Supernova. This 
formula assumes a power-law distribution of cosmic-ray energies with a typical differential 
photon index of 2.1. Calculations of the photon index give a range of values around the 
index value 2.0 for a uniform circumstellar medium (Berezhko & Volk 1997) and a rather 
more extended range of photon index values for Supernova explosions into the wind bubble 
of a very massive progenitor star (Berezhko & Volk 2002). 

A typical value for the fraction of kinetic energy converted into cosmic rays is G i^sN = 
1.0 X 10™ erg (Malkov & Drury 2001; Volk, Berezhko & Ksenofontov 2003). For an SNR at a 
distance oid — 1 kpc expanding into an environment with a typical density of the interstellar 
medium of n = 1 hydrogen atom per cm^, a VHE 7-ray flux of 20% of that from the Crab 
Nebula is expected. The luminosity Lq of such a source is 7.2 x 10^^ erg s~^ in the energy 
range between 0.2 TeV and 10 TeV. One would therefore expect that with the sensitivity 
level of the H.E.S.S. Galactic Plane survey (down to ~ 2% of the Crab flux) SNRs should 
be detectable out to a distance of 3 kpc (to be compared to 8.5 kpc to the centre of the 
Galaxy). For larger values of n, SNRs should be visible within most of the volume of the 
Galactic Disc. 

Another potential source of VHE 7-rays are pulsars and PWN. Pulsars, rapidly ro- 
tating neutron stars left over after a supernova explosion, are highly magnetised and act 
like unipolar inductors accelerating particles. The pulsar wind interacts with the ambient 
medium, generating a shock region where particles are accelerated. Such objects may there- 
fore exhibit a pulsed component of radiation, from the immediate vicinity of the pulsar, 
together with an unpulsed component from the shock region and beyond. The Crab Nebula 
is the best-studied object of this class of particle accelerators. 

Most of these potential Galactic 7-ray sources such as SNRs and pulsars are associated 
with the formation of massive stars and therefore cluster along the Galactic Plane, especially 
concentrated in the direction towards the centre of our Galaxy. 91 SNRs and 389 pulsars 
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are catalogued within the inner 60° in Galactic longitude and 6° in Galactic latitude (Green 
2004; Manchester et al. 2005). Despite the large number of known SNRs and pulsars in our 
Galaxy, the actual number of detected VHE 7-ray sources was until now limited to only a 
handful of objects (for a review, see e.g. Weekes (2000)). 

Particle acceleration to multi-TeV energies may also be powered by OB stellar associ- 
ations, via the superbubble mechanism (see Parizot ct al. (2004) for a summary), and also, 
possibly more directly via the stellar winds of member stars and subsequent interaction (see 
e.g. Montmerle (1979); Volk, & Forman (1982); White (1985); Benagha et al. (2001)). There 
are presently two TeV sources that may be linked to such OB associations. TcV J2032+4130, 
discovered by the HEGRA collaboration (Aharonian et al. 2002a), and HESS J1303-631, a 
recent discovery by H.E.S.S. (Aharonian et al. 2005b) are in the vicinity of the OB associa- 
tions Cygnus 0B2 and Centaurus OBI, respectively. Both TeV sources remain unidentified 
(UID) in the sense that there are as yet no clear counterparts of similar size and morphology 
in other wavebands. These sources might represent a new class of objects distributed along 
the Galactic Plane. 

At lower energies, observations by the space based instrument EGRET in the GeV 
energy range (Hartman et al. 1999) also revealed a population of sources clustering along 
the Galactic Plane. The limited angular resolution of the EGRET instrument renders an 
identification of these sources in other wavebands very difficult and most of these objects 
remain as yet unidentified. 

For all these reasons, a survey of the inner part of the Galaxy provides an efficient way 
to search for 7-ray emission in order to investigate properties of source classes and search 
for yet unknown types of Galactic VHE 7-ray emitters. 

At TeV energies the Milagro water-Cherenkov detector (Atkins et al. 2004) and the Ti- 
bet air-shower array (Amenomori et al. 2002) have been used to perform large-scale surveys. 
While these instruments have the advantage of a very wide field of view (~ 1 sterad), the 
sensitivity obtained by these surveys is rather limited, reaching a flux limit comparable to 
the flux of the Crab Nebula, ~ 3 x 10~^° cm~^ s"^ (for E > 200 GeV), in one year of obser- 
vations. Both surveys covered ~ 27r sterad of the northern sky and revealed no previously 
unknown significant 7-ray sources. However, some evidence for 7-ray emission was found 
by the Milagro Collaboration from a region in the Cygnus constellation and another region 
close to the Crab Nebula (Smith et al. 2005). Recently the detection of diffuse emission from 
the Galactic Plane with Milagro has been reported (Fleysher et al. 2004). 

The HEGRA instrument was the first array of imaging air Cherenkov telescopes to be 
used to survey a part of the Galactic Plane (Aharonian et al. 2002b). The range of Galactic 
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longitudes (l) —2° < I < 85° was observed; due to the location of HEGRA in the northern 
hemisphere and the resulting large zenith angles for observations of the centre of the Galactic 
Plane, the sensitivity was reduced for the central part of the Galaxy and the energy threshold 
of the observations ranged from 500 GcV to 7 TeV. No sources of VHE 7-rays were found 
in this survey, and upper limits between 15% of the Crab flux for Galactic longitudes / > 30° 
and more than 30% of the Crab flux in the inner part of the Milky Way were derived. 

Until the completion of the High Energy Stereoscopic System (H.E.S.S. ) (Hinton 2004) 
in early 2004, no VHE 7-ray survey of comparable sensitivity of the southern sky, or of the 
central region of the Galaxy has been performed. In this paper we describe a survey of this 
region with H.E.S.S. between May and July 2004 at a flux sensitivity down to 2% of the Crab 
flux. The initial results from the survey including the description of the discovery of eight 
new sources have been summarised in Aharonian et al. (2005a). This work presents detailed 
results on these eight sources. In addition we report the detection of six additional weaker 
sources of VHE 7-rays extending to lower flux levels than our previous analysis. Spectral 
information and a more detailed study on the morphologies for all sources is provided. The 
paper is structured as follows: after a description of the H.E.S.S. instrument (section 2) and 
the Galactic survey dataset used in this paper (section 3), the data analysis technique is 
presented in section 4. Section 5 gives a detailed description of the individual sources along 
with possible counterparts in other wavebands, followed by a discussion of the results in 
section 6 and some conclusions in section 7. 



2. The H.E.S.S. Instrument 

H.E.S.S. is a system of four imaging atmospheric Cherenkov telescopes, located in the 
Khomas Highlands of Namibia at a height of 1800 m above sea level. H.E.S.S. has unprece- 
dented sensitivity in the energy range from 100 GeV to several tens of TeV. The imaging 
Cherenkov technique exploits the cascade of secondary particles produced in interactions of 
primary 7-rays with molecules in the Earth's upper atmosphere. These secondary particles 
emit Cherenkov hght, which can be imaged by the Cherenkov telescopes. The energy and 
direction of the primary 7-rays can be reconstructed based on the intensity and shape of 
the detected Cherenkov light distribution. 

Multiple views of the same particle shower with several telescopes improve background 
rejection capability and the angular and energy resolution of the instrument, in comparison 
to a single telescope. Therefore, H.E.S.S. observations are performed in a stereoscopic mode, 
whereby at least two of the four telescopes must be triggered by a shower for an event to 
be recorded (Funk et al. 2004). The telescopes each have a mirror area of 107 m^ (Bernlohr 
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et al. 2003; Cornils et al. 2003). The Cherenkov light is focussed onto fast 960 pixel photo- 
multiplier cameras, which record the Cherenkov light images (Vincent et al. 2003). The 5° 
diameter field of view of the H.E.S.S. cameras is the largest of all imaging Cherenkov detec- 
tors and is especially suited for the study of extended sources and the search for unknown 
sources in surveys. This was part of the original design concept for the H.E.S.S. telescope 
system (Aharonian et al. 1997a,b). 

The H.E.S.S. instrument has an energy threshold of ~ 100 GcV at zenith before selec- 
tion cuts, an angular resolution of ~ 0.1° for individual 7-rays and a point source sensitivity 
of < 2.0 X 10~^^ cm~^ s~^ above 1 TcV (1% of the flux from the Crab Nebula) for a 5a de- 
tection in a 25 hour observation. The sensitivity of H.E.S.S. and the capability to map large 
extended sources has been demonstrated by the resolved images of the SNRs RX J1713.7- 
3946 (Aharonian et al. 2004a) and RX J0852.0-4622 (Aharonian et al. 2005c). The serendip- 
itous detections of the unidentified source HESS J1303-631 in the field of view of the binary 
pulsar PSRB1259-63 and of the composite supernova remnant G 0.9+0.1 (Aharonian et al. 
2005d) in the field of view of the Galactic Centre source HESS J1745-290 (Aharonian et al. 
2004b) demonstrate the advantage of a large flcld of view. This large fleld of view implies 
that some sources are observed up to 2 degrees off the optical axis of the telescope system. 

To illustrate the off-axis performance within the field of view of H.E.S.S., Figure 1 
summarises some basic parameters as a function of the offset. The left panel shows the 
7-ray efficiency relative to on-axis pointing, derived from Monte-Carlo simulations for two 
different zenith angles along with measurements of the Crab Nebula. The plot shows a 
weak dependence of the 7-ray efficiency on zenith angle and reasonable agreement between 
Monte-Carlo simulations and data. The right panel of Figure 1 shows the angular resolution 
of H.E.S.S. (specified by the 68% containment radius of the point spread function (PSF)) 
versus offset in the camera for observations at two different zenith angles. It is evident that 
the angular resolution is better than 0.1° and is reasonably constant in the inner 1.5° of the 
field of view. 



3. The H.E.S.S. Galactic Plane survey 

The survey was conducted in the Galactic longitude band ±30° around 1 — 0°. Runs 
of 28 minutes duration were taken at pointing positions with a spacing of 0.7° in longitude 
in three strips in Galactic latitude: b = —l°,b = 0° and b — +1°. For the data analysis 
runs were selected based on weather and hardware conditions. 95 hours of quality-selected 
data were taken in the scan mode between May and July 2004. The observation schedule 
was optimised so that individual pointings were conducted near culmination, resulting in 
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Fig. 1. — The off-axis performance of the H.E.S.S. instrument. The left figure shows the 
7-ray efficiency relative to that for an on-axis source for Monte-Carlo simulations (dashed 
line for 20° zenith angle and solid line for 45° zenith angle) and for H.E.S.S. observations of 
the Crab Nebula at a mean zenith angle of 47° (triangles). The Monte-Carlo simulations arc 
generated assuming a Crab-like source spectrum of the form dN/dE oc with F = 2.5. 
The right figure shows the simulated angular resolution (68% containment radius) as a 
function of the offset for observations at zenith angles of 20° (circles) and 50° (triangles), 
respectively. In each plot the arrow denotes the offset of 2° up to which 7-rays are included 
in this analysis. 
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the smallest zenith angles accessible. Promising source candidates were re-observed between 
July and October 2004, yielding a further 30 good-quality hours of data. For the analysis 
described in the following, additional pointed observations of SNR RX J 171 3. 7-3946 and 
of the Galactic Centre region are included. The total data set presented here amounts to 
230 hours after quality selection, the mean zenith angle of all observations is 26°. Figure 2 
shows a map of exposure for this data set, expressed in terms of an equivalent number of 
hours observing a source at the centre of the field of view. As can be seen from Figure 2, the 
exposure is not uniform with peak exposures in the regions around RX J1713. 7-3946 and the 
Galactic Centre. This irregular exposure results in an uneven sensitivity for the detection 
of new sources. The sensitivity of the survey in Galactic latitude (averaged over the whole 
longitude range) is shown in Figure 3. For —1.5° < 6 < 1.5° the sensitivity is rather flat at 
a level of 2% of the flux from the Crab Nebula, and deteriorates rapidly at > 2° from the 
Galactic Plane. The flux sensitivity above 200 GeV is ~ 50% worse for a source with photon 
index F = 3.0 compared to one with F = 2.0. 



4. Data Analysis 

Following calibration of the data (Aharonian et al. 2004c), the standard H.E.S.S. event 
reconstruction scheme was applied to the dataset. After tail-cut image cleaning and a Hillas 
parametrisation algorithm, the direction of the shower axis was reconstructed by intersection 
of the major axes of the Hillas ellipses. H.E.S.S. standard cuts (Aharonian et al. 2005e) on 
the width and length of these ellipses (scaled to the width and length expected from Monte- 
Carlo simulations) were used to suppress the hadronic background. An additional cut on the 
image size > 200 photo-electrons (p.e.) yields optimum sensitivity for weak sources provided 




Fig. 2. — The effective exposure of the H.E.S.S. Galactic Plane survey, expressed as equiva- 
lent number of hours spent observing a source at the centre of the field of view. The deeper 
observations of the RXJ1713. 7-3946 and the Galactic Centre are apparent. 
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Fig. 3. — Sensitivity in Galactic latitude of the H.E.S.S. Galactic Plane scan in units of the 
flux from the Crab Nebula (~ 3 x 10~^° cm~^ s~^, E > 200 GeV) and assuming a power- law- 
spectrum with photon index F = 2.5. The sensitivity used here is defined as the minimum 
flux of a point-source that would be detectable in this survey by H.E.S.S. at the G.Scr level 
(corresponding to Act after accounting for all trials), given the cuts used to detect slightly 
extended sources (in particular a O"^ = 0.05 deg^ cut). The solid curve shows the sensitivity 
averaged over the whole longitude range of the survey, including runs on the Galactic Centre 
region and on the SNR RXJ1713. 7-3946. The dashed curve shows the sensitivity in the 
range of longitude 350° < I < 357°, a region where no re-observations were conducted. 
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that source spectra are not significantly steeper than that of the Crab Nebula. They result 
in an angular resolution of 0.08° (68% containment radius) and an energy threshold of 
250 GeV. To reduce systematic effects which affect images close to the edge of the camera, 
only events reconstructed within a maximum distance of 2° from the pointing direction of 
the system were used for this analysis. The survey region was divided (a priori) into a grid 
of trial source positions with grid points spaced by 0.04°, a value well below the angular 
resolution of the instrument. For each grid point in the scan region, the on-source signal was 
calculated by summing up all events within a circle of radius 9 — 0.22°, (or 9'^ — 0.05 deg^), 
where 9 denotes the angular distance between the reconstructed event direction and the 
grid point. This cut is appropriate for slightly extended sources with respect to the angular 
resolution of the instrument and was chosen a-priori, since it is optimal for the detection 
of 7-ray sources with the average observed size of Galactic SNRs (Green 2004). A search 
for very extended sources (relative to the angular resolution) using a larger 9^ cut yields 
similar results. A similar search for point-like sources with a 9'^ = 0.01 deg^ cut revealed a 
subset of the sources presented here and, additionally, the point-like source HESS J 1826-148, 
which is associated with the microquasar LS5039. The discovery of this source is reported 
elsewhere (Aharonian et al. 2005f). 

The background for each grid point was estimated using two different techniques. The 
standard H.E.S.S. background estimation method counts background events in a ring around 
each grid point with a mean radius of 0.5° and an area 7 times that of the on-region to es- 
timate the background level. The standard ring radius of 0.5° was increased to 0.6° (again 
a priori) to improve sensitivity to extended sources. Using appropriate weight factors, this 
background determination takes into account the difference in system acceptance and ex- 
posure for the on- and the off-regions (Aharonian et al. 2005e). An alternative technique 
uses the template approach, where the background level is estimated using events that are 
normally rejected by the cosmic-ray background rejection criterion based on the scaled image 
width and length, but with reconstructed directions falling within the on-source region (Row- 
ell 2003) . The normahsation in this case is determined by the overall ratio of events in both 
cut regimes, taking into account the ratio of relative acceptances in the two regimes. Figure 
4 shows the correlation of the background level estimated for each grid point from the ring 
background and from the template background method. The slope of the major axis of the 
distribution is 1.007. This shows that the background estimates are consistent within 1%. 
The template background seems to slightly underestimate the background in comparison to 
the ring background. The ring background technique should therefore give a more conserva- 
tive estimation of the significances of the new sources. In the following the ring background 
estimation method will be used as the standard method, since it is a robust method in which 
linear gradients in the system acceptance cancel out. From the on-counts, the background 
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Fig. 4. — Background level estimated for each grid point, calculated using a ring background 
technique, plotted against the background level estimated from a template method (Rowell 
2003). The black dashed line shows the expected one-to-one correlation. 
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level, and the relative background acceptance in the system, a statistical significance can be 
determined for each grid point in the map, testing the existence of a 7-ray source at the 
position of that grid point (using the approach of Li & Ma (1983)). A number of significant 
sources are detected and Table 1 gives the on-counts, the background level as derived from 
the ring background determination method (after accounting for the difference in acceptance 
and in the area between the on-source and the off-source region) and the resulting number of 
7-ray excess events. Since the map contains a large number of test positions, this significance 
has to be corrected according to a "trial factor" . This trial factor accounts for the increased 
probability of finding a fake signal with an increased number of test positions for which a 
significance is calculated. Using a Monte-Carlo simulation of the technique apphed to the 
survey data (including a realistic exposure map and randomly generating the events for each 
run from the acceptance function of that run), it has been found that this trial factor is 
somewhat lower than the number of points for which a significance is calculated, i.e. the 
number of grid points in the map (n = 220000), since the excess counts for adjacent grid 
points are correlated, given that the integration radius is larger than the grid spacing. For 
uncorrelated grid points one would have exactly the factor n. For the calculation of post- 
trials significances we conservatively assume that the number of trials is equal to the number 
of bins n. Figure 5 shows the distribution of significances for all grid points in the map 
(solid black). The large positive tail to the distribution can be attributed to the existence 
of 7-r ay-sources. The much smaller tail of negative significances is due to the existence of 
7-ray sources in the regions used for background estimation. This effect has been minimised 
in a second iteration in which regions of significant 7-ray excess are excluded as background 
regions. Also shown in dashed blue is the distribution of significances for a region in the 
map without significant sources {I — ± 20°, 0.5° < b < 1.5°). This distribution has been 
scaled to have the same peak height as the distribution of significances for all grid points 
(black). A Gaussian fit to this distribution yields a mean of —0.050 ± 0.003 with an RMS of 
1.027 ± 0.003, in tolerable agreement with the expected normal Gaussian distribution given 
the possibility of undetected 7-ray-emission in the on-source or off-source regions. The dis- 
tribution of significances for a Monte-Carlo simulation of the survey without 7-ray sources, 
using the same analysis procedure as for data, is shown with red points. This distribution 
has also been scaled to have the same peak height as the distribution of significances for 
all grid points (black) and is consistent with a normal Gaussian, demonstrating that the 
analysis technique behaves as expected. A post-trial significance S2 (the effective detec- 
tion significance) was calculated from the pre-trials significance Si by converting Si into a 
probability, recalculating the probability of the occurrence of one or more such events given 
the number of trials and converting this probability back into a significance. Grid points 
for which the post-trial significance is above 4.0(7 (corresponding to 6.3a before trials) are 
considered as likely sources of VHE 7-rays and are included in the following discussion. The 
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Tablc 1: Tabic of on-source counts, scaled background level (ax A^^ofr according to the notation 
of Li & Ma (1983)) and 7-ray excess events for the new sources of VHE 7-rays detected in 
the survey of the Galactic Plane. The on-counts are calculated by counting events within 
a circle of radius 9 — 0.22° centered on the position of the source. The background counts 
are derived by counting events in a ring around the source position with average radius 0.6° 
multiplied with the normalisation factor a ~ 1/7 that takes into account the difference in 
the acceptance and the area between the on-source and the background region. 



Name 


On-counts 


Scaled Background 


Excess 


HESS J1614-518 


517 


283 


234 


HESS J1616-508 


766 


356 


410 


HESS J1632-478 


232 


139 


93 


HESS J1634-472 


275 


175 


100 


HESS J1640-465 


611 


298 


313 


HESS Jl 702-420 


311 


210 


101 


HESS J1708-410 


955 


728 


227 


HESS J1713-381 


917 


728 


189 


HESS J1745-303 


1138 


927 


211 


HESS J1804-216 


887 


508 


379 


HESS J1813-178 


717 


374 


343 


HESS J1825-137 


406 


232 


174 


HESS J1834-087 


478 


287 


191 


HESS J1837-069 


622 


328 


294 
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Significance 

Fig. 5. — Distribution of significances in tlie survey region (solid black). For comparison, 
the distribution of significances in a region without 7-ray sources is shown in dashed blue. 
The distribution of significances for a Monte- Carlo simulation of the survey without 7-ray 
sources is shown as red points and follows a Gaussian distribution with zero mean and unit 
width. A Gaussian fit to the dashed blue distribution yields a mean of —0.050 ± 0.003 with 
an RMS of 1.027 ± 0.001, indicating a small systematic error in the determination of the 
significances at the level of a few percent. The dashed blue curve and the red data points 
have been scaled to the peak height of the black distribution. 
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probability of a single fake signal at this level is 3 x 10""^ in the whole survey. Figure 6 shows 
the map of significances before accounting for the trial factor with the background estimated 
from the ring technique. The strong 7-ray sources RXJ1713. 7-3946 and HESS J1745-290 
clearly stand out, partly because of the high exposure in this regions. Also visible at a much 
fainter level is the recently published new 7-ray source G 0.9+0.1 (Aharonian ct al. 2005d), 
which was serendipitously detected in the H.E.S.S. observations of the Galactic Centre re- 
gion. In addition to these three known objects, fourteen new sources of VHE 7-rays are 
detected at a significance level above 4(7 (post-trials), the eight most significant of which 
have been discussed in a previous publication (Aharonian et al. 2005a). Table 2 shows the 
pre-trials significance for the peak positions in the map for the new sources in column Si 
and the significance after accounting for the number of trials in column 5*2. Column 5*3 
shows the post-trial significance applying a cut appropriate for point sources (^^point = 0.1°) 
instead of the extended source cut 6 = 0.22°. Only two of the new sources (HESS J1640-465 
and HESS J1813-178) show a non-marginal increase in significance using the point-source 
cut . The extension of the new sources will be discussed in the next section. It should be 
noted, that the five sources with detection significances yS'2 below 5a should be considered as 
source candidates, rather than firmly establish sources. Reobservation of these sources will 
be necessary to confirm these marginal detections. 

5. Chciracteristics of the new sources 

This section describes the main results and the analysis techniques that were applied to 
provide information about position, extension and spectra of the new VHE 7-ray sources. 

5.1. Position and Morphology 

A detailed spectral and positional analysis has been performed on the fourteen sources 
listed in Table 2. For this purpose, the distribution of 7-ray candidates was determined for 
fine angular bins of width 0.005°, to allow for a more accurate determination of the source 
position. For the position and approximate size determination of each source, the spatial 
distribution of events in the (uncorrelated) bins was fit to a model of a 2-D Gaussian 7-ray 
brightness profile of the form p oc exp{—9'^/2a'^^^^^), convolved with the PSF of the instru- 
ment, and an added background parametrisation. To properly handle the small number 
of events in the bins, a likelihood fit appropriate for Poisson statistics was applied in this 
procedure. To test for possible elongation of each emission region, an elongated Gaussian 
with independent crgource and cr'sonrce ^Sich dimension and a free orientation angle uj (mea- 
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Galactic Longitude (°) 



Fig. 6. — Significance Map of tlie H.E.S.S. Galactic Plane survey in 2004, including data 
from re-observations of source candidates detected in the original scan and observations of 
the known 7-ray emitter RX J 171 3. 7-3946 and the Galactic Centre region. The typical 
energy threshold for this map is 250 GeV. The on-source counts for each grid point are 
integrated in a circle of radius 9 = 0.22°. The background for each grid point has been 
derived using a ring of mean radius 0.6° and an area 7 times that of the on-source circle. 
The labels indicate the 7-ray sources described in this work, along with the known 7-ray 
sources RX J 171 3. 7-3946 (HESS J1713-397), G 0.9+0.1 (HESS J1747-281) and the Galactic 
Centre (HESS J1745-290). The numbers in the map give the post-trials significances of the 
7-ray sources. The significance scale is truncated at 18a, the signals from the Galactic 
Centre and RX J1713. 7-3946 exceed this level. 
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Table 2: Table of significances for the new sources of VHE 7-rays detected in the survey 
of the Galactic Plane. Source positions are derived from a fit to the 7-ray counts in fine 
angular bins, and do not rely on the coarser spaced search grid. Si stands for the pre-trials 
significance (using a 9"^ cut of 0.05 deg^) for the peak positions in the map. S2 shows the 
post-trial significance for this position. Column S3 shows the post-trial significance for this 
position for a point-source cut (^point = 0.1°). The sources that are marked with a ^ have 
been previously reported in (Aharonian et al. 2005a). Galactic longitude I and latitude b as 
well as RA/Dec coordinates are given in degrees, the live time is given in hours. 



Name 








Best Fit Position 


Sig 


;nificance 


Live time 








in 


bn 


RA 


n 


Dec(°) 


Si 




^3 




HESS 


J1614- 


-518^ 


331.52 


-0.58 


243 


58 


-51.82 


11.2 


10.1 


7.5 


9.8 


HESS 


J1616- 


-508t 


332.39 


-0.14 


244 


10 


-50.90 


17.1 


16.3 


9.7 


10.2 


HESS 


J1632- 


-478 


336.38 


0.19 


248 


04 


-47.82 


6.6 


4.5 




4.5 


HESS 


J1634- 


-472 


337.11 


0.22 


248 


74 


-47.27 


6.4 


4.2 


1.8 


6.6 


HESS 


J1640- 


-465^ 


338.32 


-0.02 


250 


18 


-46.53 


14.5 


13.6 


14.5 


14.3 


HESS 


J1702- 


-420 


344.26 


-0.22 


255 


69 


-42.07 


6.3 


4.0 


4.1 


5.7 


HESS 


J1708- 


-410 


345.67 


-0.44 


257 


06 


-41.08 


8.4 


6.8 


6.9 


37.4 


HESS 


J1713- 


-381 


348.65 


0.38 


258 


49 


-38.20 


6.3 


4.0 


1.4 


37.3 


HESS 


J1745- 


303 


358.71 


-0.64 


266 


26 


-30.37 


6.3 


4.0 


1.8 


35.3 


HESS 


J1804- 


-216t 


8.40 


-0.03 


271 


13 


-21.70 


13.9 


13.0 


7.7 


15.7 


HESS 


J1813- 


-178t 


12.81 


-0.03 


273 


40 


-17.84 


14.4 


13.5 


15.3 


9.7 


HESS 


J1825- 


-137t 


17.82 


-0.74 


276 


51 


-13.76 


9.5 


8.1 


3.7 


8.4 


HESS 


J1834- 


-087^ 


23.24 


-0.32 


278 


69 


-8.76 


9.5 


8.1 


5.7 


7.3 


HESS 


J1837- 


-069t 


25.18 


-0.11 


279 


41 


-6.95 


13.2 


12.2 


9.9 


7.6 
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sured counter-clockwise from the positive Galactic latitude axis) was also fit. The source 
is flagged as elongated if the error range in best fit minor axis cr'^^^^^e is outside that of the 
best fit major axis cXsource- This fitting procedure has been tested using Monte-Carlo sim- 
ulations under realistic conditions. The re-analysis of the source positions presented here 
in comparison to (Aharonian et al. 2005a) yields a maximum shift of 2-3 arc minutes for 
the most extended of these sources, reflecting the more detailed study presented here. The 
new positions are consistent within errors with the previously published positions. Also the 
previously reported extensions are confirmed by the more detailed analysis presented here, 
with the exception of HESS J1834-087, where the extension was overestimated due to the 
non-Gaussian shape of the source and is now determined to be 0.09° rather than 0.20°. 

The system PSF (see Figure 1) has been derived from Monte Carlo simulations as a 
function of the zenith angle of observations and of the offset of a source in the field of view. 
This model PSF is consistent with the excesses observed for the point-hke Crab Nebula 
and the active galaxy PKS 2155-304 (Aharonian et al. 2005e). Table 3 shows the best 
fit position, the best fit Gaussian equivalent sizes o-gourcc and (Tso^j.^.^, and the orientation 
00 (in case of elongation) of the fourteen sources. To quantify the quality of the fit, the 
last column of Table 3 shows the normalised between the model source shape and the 
measured radial profile fo each source. Figure 7 shows the distribution of sizes for the new 
sources. For the elongated sources, the arithmetic mean of cXsource and o-'^omce been filled 
into the histogram. All but two of the sources appear significantly extended above a level 
of three standard deviations. Four out of fourteen appear significantly elongated, whereas 
the remaining ten new sources are compatible with a radially symmetric shape. The radial 
excess distributions of all sources are given individually in section 5.3. 



5.2. Spectra 

Using the best fit position, a spectral analysis was performed, again assuming a Gaus- 
sian emission region and using a cut of 9^ < ^point ("^source + '^psf)/'^psf' approximating the 
H.E.S.S. PSF as a 2-D Gaussian with ctpsf = 0.08° and ^point = 0.1° being the optimal cut 
for point sources. Such a cut provides optimum significance for a weak Gaussian shaped 
source on a uniform background, and hence the best possible spectral determination. For 
the elongated sources, the larger of the two values (Xsource and cr'sonicc used. The en- 
ergy dependent efficiency of this 6 cut is estimated by convolving the H.E.S.S. PSF with a 
Gaussian of size cXsource- This introduces an additional systematic error in the absolute flux 
determination since it implicitly assumes that the emission region of the source follows a 
Gaussian distribution and that the morphology of the source does not show any energy de- 
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Table 3: Position and sizes of the new sources of VHE 7-rays detected in the survey of the 
Galactic Plane. Galactic longitude I and latitude b, the RMS source size (Jsource and crgguj.^e 
and the orientation angle cu counterclock-wise with respect to the positive Galactic latitude 
axis are given in degrees. Xnorm denotes the normalised x^, with 49 degrees of freedom. 



Name 




Best Fit Position 




''"source 


UJ 


A norm 






in 


bn 


n 


n 


n 




J1614- 


-518 


331.52±0.03 


-0.58±0.02 


0.23±0.02 


0.15±0.02 


49±10 


1.44 


J1616- 


-508 


332.391±0.014 


-0.138±0.013 


0.136±0.008 






0.82 


J1632- 


-478 


336.38±0.04 


0.19±0.03 


0.21±0.05 


0.06±0.04 


21±13 


1.05 


J1634- 


-472 


337.11±0.05 


0.22±0.04 


0.11±0.03 






0.72 


J1640- 


-465 


338.316±0.007 


-0.021±0.007 


0.045±0.009 






0.91 


J1702- 


-420 


344.26±0.04 


-0.22±0.03 


0.08±0.04 






1.38 


J1708- 


-410 


345.672±0.013 


-0.438±0.010 


0.054±0.012 






0.86 


J1713- 


-381 


348.65±0.03 


0.38±0.03 


0.06±0.04 






1.04 


J1745- 


-303 


358.71±0.04 


-0.64±0.05 


0.21±0.06 


0.09±0.04 


54±7 


1.17 


J1804- 


-216 


8.401 ±0.016 


-0.033±0.018 


0.200±0.010 






1.34 


J1813- 


-178 


12.813±0.005 


-0.0342±0.005 


0.036±0.006 






0.94 


J1825- 


-137 


17.82±0.03 


-0.74±0.03 


0.16±0.02 






1.03 


J1834- 


-087 


23.24±0.02 


-0.32±0.02 


0.09±0.02 






1.16 


J1837- 


-069 


25.185±0.012 


-0.106±0.016 


0.12±0.02 


0.05±0.02 


149±10 


1.52 
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Fig. 7. — Distribution of the RMS angular size (Xsourcc of the fourteen VHE 7-ray sources. 
For the elongated sources, the arithmetic mean of agource and o-'g^urce used. The mean of 
this distribution is 0.1° with an RMS of 0.05°. 
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pendency. The typical overall systematic error in the determination of the integral flux level 
is estimated to be ~ 30%. The standard H.E.S.S. spectral analysis technique is described 
in (Aharonian et al. 2005e). To avoid systematic errors in the spectral determination arising 
from different camera acceptances across the H.E.S.S. field of view, the spectral background 
determination employs regions of the same size and shape as the on region, but displaced 
on a ring in the field of view around the camera centre. The ring was chosen such that the 
off regions have the same angular distance to the centre as the on region. This approach 
ensures that background events are taken at the same zenith and offset angles. 

To keep as many low-energy 7-rays as possible, a looser size cut of 80 photo-electrons 
is used in the spectral analysis. This cut allows to reach a lower energy threshold (150 GeV) 
at the expense of reduced sensitivity for hard-spectrum sources and a 30% poorer angular 
resolution. Flux points below a significance level of 2 standard deviations were discarded 
in the spectral analysis and the binning was doubled at high energies for weak sources. 
Tabic 4 lists the resulting photon index F and flux above 200 GeV for a power-law fit 
dN/dE = FqE^^. The typical systematic error on these photon indices is estimated to be 
~ 0.2. It should be noted, that the fluxes reported here are in general different from the 
fluxes reported in Aharonian et al. (2005a) (especially for the larger sources) since in the 
present work the flux is determined from a larger region, adapted to encompass a larger 
fraction of the emission region, whereas in Aharonian et al. (2005a), a fixed small 9"^ — 0.02 
deg^ cut was used for all sources. Figure 8 shows the distribution of photon indices for the 
fourteen new sources. The distribution exhibits a mean of 2.32, marginally consistent with 
predictions by models of shock wave acceleration in SNRs (see e.g. Berezhko & Volk (1997, 
2002)) and as expected for the source spectrum of Galactic cosmic rays (see e.g Jones et al. 
(2001)). The distribution has an RMS of 0.2, comparable with the statistical errors on the 
measurements. 



5.3. Individual 7-ray sources 

This section describes the fourteen VHE 7-ray sources detected in the Galactic Plane 
survey individually along with a brief description of possible counterparts at other wave- 
lengths. Excess maps are shown that have been smoothed with a two-dimensional Gaussian 
to reduce the effect of statistical fluctuations. The resulting count map is in units of inte- 
grated excess counts within the smoothing radius. The RMS radius of the Gaussian used for 
the smoothing is adapted for each source according to the available photon statistics. This 
radius is given in the flgure captions, and is shown as a black dotted line in the inset in the 
lower right part of each excess map along with the PSF of the instrument, smoothed in the 
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Table 4: Results of a power- law fit dN/dE — FqE~^ for the fourteen sources of VHE 7-rays. 
The d cut used for the spectral analysis is given in degrees. The errors given on the photon 
index and flux are statistical only. The typical systematic errors are ~ 0.2 in the photon 
index and ~ 30% in the absolute flux. To derive a combined error on the photon index, the 
systematic and statistical errors should be added in quadrature. 



Name 




r 


Flux > 200 GeV 

{!{)-'- cm-2 s-i) 


9 cut 

(°) 


XVNDF 


HESSJ1614- 


-518 


2.46±0.20 


57.8±7.7 


0.40 


0.6/3 


HESS J1616- 


-508 


2.35±0.06 


43.3±2.0 


0.26 


5.5/5 


HESS J1632- 


-478 


2.12±0.20 


28.7±5.3 


0.36 


7.9/2 


HESS J1634- 


-472 


2.38±0.27 


13.4±2.6 


0.23 


5.8/5 


HESS J1640- 


-465 


2.42±0.15 


20.9±2.2 


0.16 


3.7/4 


HESSJ1702- 


-420 


2.31±0.15 


15.9±1.8 


0.20 


7.8/6 


HESSJ1708- 


-410 


2.34±0.11 


8.8±0.7 


0.17 


6.3/3 


HESSJ1713- 


-381 


2.27±0.48 


4.2±1.5 


0.17 


1.8/2 


HESSJ1745- 


-303 


1.82±0.29 


11.2±4.0 


0.36 


4.4/3 


HESS J1804- 


-216 


2.72±0.06 


53.2±2.0 


0.36 


6.9/11 


HESS J1813- 


-178 


2.09±0.08 


14.2±1.1 


0.15 


5.5/11 


HESS J1825- 


-137 


2.46±0.08 


39.4±2.2 


0.30 


6.2/8 


HESS J1834- 


-087 


2.45±0.16 


18.7±2.0 


0.20 


3.5/4 


HESSJ1837- 


-069 


2.27±0.06 


30.4±1.6 


0.23 


12.6/10 
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same way as the excess map. Note that the maps have different scales, therefore the white 
hne in the upper right corner indicates an 0.2° scale. Further it should be noted that the 
colour scale covers different intervals for different maps. 

Also shown are possible counterparts to the VHE 7-ray emission: SNR from Green 
(2004) (with white circles illustrating their nominal radii), pulsars from the ATNF pulsar 

catalogTic (Manchester et al. 2005) (as white triangles), hard X-ray soiuxcs detected by the 
International Gamma- Ray Laboratory (INTEGRAL) satellite (as white stars), and uniden- 
tified EGRET sources from the third EGRET catalogue (Hartman et al. 1999) (showing the 
95% error contour for the position of the source in dashed white). In some maps contours of 
radio and/or X-ray emission are also shown and discussed individually. 20 cm Very large ar- 
ray (VLA) radio contours are taken from White, Becker & Helfand (2005). The x-ray ASCA 
data show smoothed count maps obtained with the GIS instrument in the 0.7-10 keV band. 
The best fit source position of the H.E.S.S. source is shown with error bars in black, the best 
fit (Gaussian) source size is shown as a black dashed circle or ellipse. The radial distribution 
of 7-ray excess with respect to the best-fit source position is shown for each source binned 
appropriately for the available photon statistics; in each plot the number of reconstructed 
events in bins of the squared angular distance {6^) is shown along with the best-fit radially 
symmetric model (solid black line) and the PSF of the instrument (dashed red line), derived 
from Monte-Carlo simulations. For all sources the reconstructed differential energy spectra 
are shown together with the best- fit power-law (solid line). 



5.3.1. HESSJI6I4-5I8 

This source (see Figure 9) is located in a region relatively devoid of counterpart candi- 
dates. The source exhibits elliptical emission with a semi major axis of 14±1 arc minutes and 
a semi minor axis of 9±1 arc minutes. HESS J1614-518 is the brightest of the new sources, 
detected with a fiux corresponding to 25% of that from the Crab Nebula above 200 GeV. 
The spectrum is well fit by a power- law of photon index 2.46±0.20, see Figure 9 (top right). 
A nearby Chandra observation (Kastner et al. 2003) covered an area more than 10 arc min- 
utes away from the best fit position of HESS J1614-518 and showed no evidence for X-ray 
emission from the direction of the H.E.S.S. source. The poor quality of the fit to the radial 
distribution of HESS J1614-518 as is shown in Figure 9 (bottom right) indicates that the 
source is not well described by a symmetrical Gaussian profile. 
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Fig. 8. — Distributions of the photon index of the new sources. The mean photon index is 
2.32 with an RMS of 0.2. 
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Fig. 9. — Smoothed excess map (left, smoothing radius 0.09°) of the region of HESS J1614- 
518. The figures on the right hand side show the energy spectrum (top) and 6"^ plot (bottom). 
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5.3.2. HESSJ1616-508 

HESS J1616-508 (see Figure 10) with a flux of 19% of tlie Crab Nebula above 200 GeV 
is also one of the brighter new sources of VHE 7-rays. Its spectrum can be fitted by a power- 
law of photon index 2.35±0.06 as shown in Figure 10 (top right). The emission region has 
an RMS size of 16 arc minutes and does not have a good positional match to any plausible 
counterpart at other wavelengths. It is located in a region close to the known SNRs G332.4- 
0.4 (ROW 103) and G332.4+0.1 (Kes32). The region has been intensively studied with 
X-ray satellites like Chandra (Vink 2004) and XMM, with the position of HESS J1616-508 
mostly outside or at the edge of the field of view where the sensitivity is strongly decreased. 
No extended X-ray source that lines up with the VHE 7-ray emission has been found. 
Another interesting object in this region is the energetic pulsar PSR J1617-5055 (Torii et al. 
1998), a young X-ray-emitting pulsar, that was proposed to emit TeV 7-rays (Hirotani 2001). 
Approximately 1.2% of its spin-down luminosity would suffice to power the observed (0.2-10 
TeV) VHE 7-ray flux, assuming the H.E.S.S. source is located at the pulsar distance. Since 
HESS J1616-508 does not fine up with this object, an asymmetric (yet undetected) PWN 
nebula powered by the pulsar would be required to explain the TeV emission. 



5.3.3. HESS J 1632-478 

The centroid position of this object (right hand source in Figure 11 (top)) and of the 
close-by HESS J1634-472 are spatially separated by less than 45 arc minutes. HESS J1632- 
478 has an elongated shape with a semi-major axis of 12±3 arc minutes and a semi-minor 
axis of 3.6±2.4 arc minutes. Its flux in VHE 7-rays above 200 GeV corresponds to 12% of 

the flux from the Crab Nebula and the spectrum can be fitted by a power-law of photon in- 
dex 2.12±0.20 (see Figure 11 (bottom right)). A positional coincidence of HESS J1632-478 
exists with the hard X-ray source IGR J16320-4751 discovered by the INTEGRAL satel- 
lite (Tomsick et al. 2003) in the energy range above 15 keV. In softer X-rays between 2 and 
10 keV, this object was also strongly detected in an XMM-Newton observation (Rodriguez 
et al. 2003) and found to be coincident with the ASGA source AX J 163 1.9-4752, detected in 
the ASGA Galactic Plane survey (Sugizaki et al. 2001). This source is beheved to belong to a 
new class of heavily absorbed high mass X-ray binary systems, with an equivalent absorption 
column density of (9.2 ±1.1) x 10^^ cm^^ in this case. The extension of HESS J 1632-478 does 
not plead in favor of an association with IGR J16320-4751. One should note that the proba- 
bility of spurious associations is increased in this region of the Galactic plane (the 3 kpc arm 
tangent region) due to the high density of absorbed INTEGRAL sources. Another potential 
counterpart is the ASCA source AX J 163252-4746. 
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Fig. 10. — Smoothed excess map (left, smoothing radius 0.06°) of the region surrounding 
HESS J1616-508 (left hand source), along with nearby pulsars and SNRs which were consid- 
ered as counterparts. HESS J1614-518 is also visible in this map (right hand source). The 
figures on the right hand side show the energy spectrum (top) and the radial distribution 
(e^) plot (bottom). 
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Fig. 11. — Smoothed excess map (top, smoothing radius 0.12°) of the region of the two close- 
by H.E.S.S. sources HESS J1634-472 (left hand source) and HESS J1632-478 (right hand 
source), along with nearby pulsars, SNRs, and INTEGRAL sources which were considered 
as counterparts. It should be noted, that the tail of HESS J1632-478 to negative Galactic 
latitude is statistically not significant. The figures in the middle row show the spectra (left: 
HESS J1634-472, right HESS J1632-478), the figures in the bottom row the corresponding 
plots. 
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5.34- HESSJ1634-472 

The size of the HESS J1634-472 (left hand source in Figure 11 (top)) emission region is 
6.6±1.8 arc minutes and its flux above 200 GeV is close to 6% of the Crab flux with a photon 
index for the energy spectrum of 2.38±0.26 (see Figure 11, centre left). HESS J1634-472 
has no direct positional counterpart, although it is interesting to note that yet another IN- 
TEGRAL source IGR J16358-4726 is located close-by (Revnivtsev et al. 2003). This source 
presumably belongs to the class of highly absorbed X-ray binary systems, like IGR J16320- 
4751 (discussed in the previous paragraph). The source was detected in the 2-10 keV X-ray 
band with the Chandra satellite (Patel et al. 2004). The column density dermined for this 
object is Nh = 3.3 x 10^^ cm~^ and the luminosity of this object in the 2-10 keV band is 
]^g36 QYg/s. Another potential counterpart of HESS J1634-472 is the recently discovered SNR 
candidate G337.2+0.1 (Combi et al. 2005), coincident with the ASCA source AXJ1635.9- 
4719. However, the distances of the two objects to HESS J1634-472 renders a correlation 
unlikely. 



5.3.5. HESSJ1640-465 

As can be seen from the radial distribution of this source (see Figure 12 (bottom right)) 
HESS J1640-465 is marginally extended with respect to the PSF of the instrument. It is 
one of the few examples of a excellent spatial correlation between the VHE 7-ray emission 
and an SNR (G338. 3-0.0), known from radio observations (Green 2004). G338. 3-0.0 was 
observed in the Molonglo Galactic Plane survey at 843 MHz (shown in green contours in 
Figure 12 (left)) and has been reported to be a broken shell SNR lying on the edge of 
a bright HII region (Whiteoak & Green 1996). The region bridging this SNR and the 
close-by SNR G 338.5+0.1 is known to contain dense HII regions. The relation of surface- 
brightness to distance yields a distance of 8.6 kpc, which would place the SNR in the 3 kpc 
arm tangent region. The ASCA source AX J164042-4632 with a photon index of 2.98!l^o 89 
and a flux level of 1.21 x 10^^^ ergs cm~^ s~^, detected during the ASCA Galactic Plane 
survey (Sugizaki et al. 2001), has been identified with G338. 3-0.0. It should be noted, 
that the position of HESS J 1640-465 is compatible with the position of the unidentified 
EGRET source 3EG J 1639-4702 (Hartman et al. 1999). The 95% positional confidence 
contour of 3EG J1639-4702 is shown in Figure 12 (left). The distance of HESS J1640-465 
to this EGRET source is 34 arc minutes. The differential energy spectrum of HESS J1640- 
465 can be well flt by a power-law with photon index 2.42±0.14 at a flux close to 9% of 
that from the Crab Nebula above 200 GeV (see Figure 12 (top right)). A simple power-law 
extrapolation of the VHE gamma-ray spectrum to 1 GeV matches the EGRET flux above 
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that energy. 

5.3.6. HESS J 1702-420 

For HESS J1702-420 (Figure 13) no significant evidence for extension of the emission 
region can be derived beyond a level of 2 standard deviations for the measured source size 
of 4.8±2.4 arc minutes. The flux above 200 GeV of this source corresponds to 7% of the 
Crab flux and the energy spectrum can be fitted by a power-law of photon index 2.31±0.15 
(see Figure 13 (top right)). The close-by pulsar PSR J1702-4128 is energetic enough to 
account for the observed VHE 7-ray emission by converting roughly 14% of its spin-down 
luminosity but would need an asymmetric PWN to be responsible for the observed VHE 
7-ray emission. No asymmetric PWN has been detected from this pulsar, however, which 
renders an association unlikely. The nearby SNR G 344.7-0.1 is too small and not positionally 
coincident with the VHE 7-ray-source. Therefore, we cannot identify a plausible counterpart 
for this object. 

5.3.7. HESSJI7O8-4IO 

This VHE 7-ray source (Figure 14) of size 3.2±0.7 arc minutes is located very close to 
SNR RXJ1713. 7-3946 and therefore its observation live time (within 2° of the pointing posi- 
tion of the instrument) amounts to 37.2 hours, well beyond typical exposures in the H.E.S.S. 
survey. The source exhibits an energy spectrum with a power-law index of 2.34±0.11, the 
flux above 200 GeV is less than 4% of the flux from the Crab Nebula (see Figure 14 (top 
right)). As can be seen from Figure 14 (left), there is no close-by plausible counterpart. 
A nearby XMM observation of the SNR G345.7-0.2 has revealed no X-ray counterpart to 
the VHE 7-ray source, but it should be noted that the offset of HESS J1708-410 from the 
centre of the XMM field of view was 15 arc minutes, positioning it close to the edge of the 
field of view. As can be seen from the radial distribution, shown in Figure 14 (bottom right), 
the source is only marginally extended with respect to the PSF of the instrument. 

5.3.8. HESSJ1713-381 

This source (Figure 15) is the second of the fourteen VHE sources that shows no 
significant evidence beyond 2 standard deviations for an extension of the emission region. 
The energy spectrum of HESS J1713-381 is consistent with a power-law with photon index 
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Fig. 12. — Smoothed excess map (left, smoothing radius 0.05°) of the region surrounding 
HESS J1640-465 along with nearby SNRs which were considered as counterparts as discussed 
in the text. The green contours indicate the radio emission detected from the region sur- 
rounding G338. 3-0.0 in the Molonglo Galactic Plane survey at 843 MHz (Green et al. 1999). 
The dashed white contours show the 95% positional error of the unidentified EGRET source 
3EG J1639-4702. The differential energy spectrum (top) and the radial distribution (6*^) plot 
(bottom) are shown on the right hand side. 
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Fig. 13. — Smoothed excess map (left, smoothing radius 0.10°) of the region surrounding 
HESS 1702-420 along with nearby pulsars and SNRs which were considered as counterparts 
as discussed in the text. The tail extending to positive Galactic longitude and latitude is not 
statistically significant. The differential energy spectrum (top) and the radial distribution 
(6*^) plot (bottom) are shown on the right hand side. 
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2.2±0.5. Its flux of 1.8% of the flux from the Crab Nebula above 200 GeV makes this object 

the weakest of the new sources. Its location is very close to the strong VHE7-ray source 
SNR RXJ1713. 7-3946 and it was therefore in the field of view of most pointed observations 
by H.E.S.S. on this object. The correspondingly deep exposure of 37 hours at offsets of less 
than 2° from the centre of the field of view made it possible to detect this object. As can be 
seen from Figure 15 (left), it coincides with part of the unusual SNR complex CTB 37. The 
radio contours of the CTB 37 complex as taken from the Molonglo Plane survey at 843 MHz 
are shown in white (Green et al. 1999). There is clear positional match between the source 
HESSJ1713-381 and the SNR G 348.7+0.3 (CTB 37 B). Also visible in Figure 15 (right) 
is a weaker marginal source consistent with the position of G 348.5+0.1 (CTB 37 A). The 
marginal detection of emission from this source is however not significant after accounting 
for all trials for the analysis of extended sources in this survey as described above. Kassim, 
Baum & Weiler (1991) describe CTB 37 A as two SNRs overlapping in projection: with the 
partial shell of G 348.7+0.1 in the west and the SNR G 348.5-0.0 to the east. Both SNRs 
in CTB 37A appear to be interacting with a system of molecular clouds with densities of 
100-1000 cm-3 (Reynoso & Mangum 2000). The SNR CTB 37 B is less weU studied but 
appears to share the dense environment of CTB 37 A. Evidence for X-ray emission from the 
CTB 37 complex can also be seen in archival ASCA data. 



5.3.9. HESSJ1745-303 

This object (Figure 16) has a direct positional coincidence with an unidentified EGRET 
source (3EG J1744-3011) (Hartman et al. 1999). Since it is within the field of view of the 
Galactic Centre, HESS J1745-303 has a large observation live time of 35.3 hours. Its flux 
level corresponds to 5% of the flux from the Crab Nebula above 200 GeV and the emission 
exhibits a hard differential ffux of index 1.8+0.3. Its emission region has an elliptical shape 
with a semi-major axis of 13+4 arc minutes and a semi-minor axis of 5.4+2.4 arc minutes as 
is shown in Figure 16 (left). The distance of HESS J1745-303 to the centroid of the EGRET 
source 3EG J1744-3011 is 10 arc minutes, which is well within the 95% uncertainty level of 
the EGRET position. A simple power-law extrapolation of the VHE 7-ray spectrum to 
energies above 1 GeV is an order of magnitude below the EGRET flux above that energy 
(F(> 100 MeV) = (64 + 7) x 10~^ cm~^ s~l). Another potential counterpart is the pulsar 
PSRB1742-30 but essentially 100% of its spin-down power would be required to account for 
the entire VHE 7-ray emission. 



5.3.10. HESSJ1804-216 



This source (Figure 17) with a size of approximately 22 arc minutes is the largest of the 
new emission regions of VHE 7-rays. With a flux of nearly 25% of the flux from the Crab 
Nebula above 200 GeV, it is also the brightest of the new H.E.S.S. sources. The photon 
index of 2.72±0.06 makes it one of the softest sources (see Figure 17 (top right)). As can 
be seen from Figure 17 (left), the emission region does not perfectly line up with any known 
object, although it could be associated with the south-western part of the shell of the SNR 
G8. 7-0.1 of radius 26 arc minutes. From CO observations, the surrounding region is known 
to be associated with molecular gas where massive star formation takes place (Ojeda-May 
et al. 2002). The white contours in Figure 17 (left) show soft X-ray emission as detected 
by the ROSAT satellite (Finley & Ogelman 1994), the black contours show the 20 cm radio 
emission as recorded with the VLA. Another plausible association of HESS J1804-216 is 
the young Vela-like pulsar PSR J1803-2137 with a spin-down age of 16000 years (Kassim 
& Weiler 1990). The required efficiency of this pulsar's spin-down luminosity to power the 
observed emission in 7-rays is only 3% and it could therefore easily account for the new 
source. 



5.3.11. HESS 31813-178 

With a moderate flux of approximately 6% of the Crab Nebula above 200 GeV, this 
source (see Figure 18) is among the more compact of the new sources with an extension of 2.2 
arc minutes. When using a large 9'^ cut of (0.4°)^, the normalisation of the power- law flt is 
25% higher, suggesting a faint non-Gaussian tail in the radial brightness proflle. The source 
exhibits a rather hard photon index of 2.09±0.08 (see Figure 18 (top right)). In the previous 
publication (Aharonian et al. 2005a) this source was still considered as unidentifled. As it is 
clear from Figure 18 (left), there is a compelling positional coincidence between this source 
and the X-ray source AX J1813-178 (Brogan et al. 2005) (white contours in flgure 18 for the 
0.7-10 keV band). This ASCA source is one of the bright sources detected in the ASCA 
Galactic Plane survey with a very hard energy spectrum. This object was also recently 
detected in the INTEGRAL Galactic plane survey in the 20-100 keV energy band (Ubertini 
et al. 2005). The region surrounding HESS J1813-178 was also observed with the VLA and 
20 cm radio contours are shown in black in Figure 18. There is a faint radio source (G 12.82- 
0.02) with a shell-like structure visible at the position of HESS J1813-178 as also recently 
reported (Brogan et al. 2005; Helfand, Becker & White 2005). This indicates a probable 
association with a shell-type SNR. Additionally HESS J1813-178 lies 10 arc minutes from 
the centre of the radio source W33. W33 extends over 15 arc minutes and has a compact 
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radio core (G 12.8-0.2) which is 1 arc minute across (Haschick Sz Ho 1983), located ~ 15 arc 
minutes away from HESS J1813-178. The W33 region shows signature of an ultra-compact 
HII region (Churchwell 1990) and contains methanol, hydroxyl and water masers and other 
tracers of recent star formation, that could act as target material for the generation of VHE 
7-rays but is not coincident with the H.E.S.S. source. 

5.3.12. HESSJ1825-137 

This source has an energy spectrum with a photon index of 2.46±0.08 and a flux of 
VHE 7-rays of approximately 17% of the flux from the Crab Nebula above 200 GeV (see 
Figure 19 (top right)). It is approximately radially symmetric and the emission region can 
be best described by a Gaussian of width 10 arc minutes. The position of HESSJ1825- 
137 is approximately 11 arc minutes to negative Galactic longitudes of the young pulsar 
PSRB1823-13 (PSR J1826-1334), which makes an association of the VHE 7-ray source 
with the pulsar possible. The pulsar PSRB1823-13 is a 101 ms evolved Vela-like pulsar 
at a distance of 3.9±0.4 kpc and at this distance the VHE luminosity of HESS J1825-137 
would only be ~2% of its spin-down power. The region surrounding the pulsar has been 
previously observed at large zenith angles by the Whipple collaboration (Hall et al. 2003), 
who observed a 3.1(T excess also to negative Galactic longitudes with respect to the pulsar 
position. X-ray observations with the XMM-Newton satellite show a 5 arc minutes diameter 
diffuse emission region extending also asymmetrically to negative Galactic longitudes from 
the pulsar (Gaensler et al. 2003). If the VHE 7-ray source is indeed associated with this 
pulsar wind nebula, the VHE emission region would extend much further to the south than 
in X-rays. PSRB1823-13 was also proposed to power the close-by unidentified EGRET 
source 3EG J1826-1302 (Nolan et al. 2003). The distance of the centroid of this EGRET 
source to HESS J1825-137 is 43 arc minutes. As illustrated in Figure 19 (left), the new VHE 
7-ray source is well compatible with the EGRET source within the 95% positional error of 
EGRET and a simple power-law extrapolation of the H.E.S.S. energy spectrum to energies 
above 1 GeV gives a similar flux as derived from the EGRET spectrum above this energy. 
An in-depth study of HESS J1825-137 and dicussion about the corresponding PWN model 
will be given in a separate paper (Aharonian et al. 2005h). 

5.3.13. HESS J1834-087 

HESS J1834-087 (Figure 20) has a size of 12 arc minutes and a 7-ray flux of 8% of the 
flux from the Crab Nebula above 200 GeV. Its energy spectrum exhibits a photon index of 
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2.45±0.16 (see Figure 20 (top right)). When using a large 9'^ cut of (0.4°)^, the normahsation 
of the power-law fit increases by 30%, suggesting a deviation of the morphology from the 
assumed radial Gaussian brightness profile. The location of HESS J1834-087 makes it an- 
other source with a compelling positional agreement with an SNR (see Figure 20 (left)). The 
positionally coincident shell-type SNR G23.3-0.3 (W41) with a diameter of 27 arc minutes 
was mapped in radio with the VLA at 330 MHz (Kassim 1992) and at 20 cm (White, Becker 
& Helfand 2005) as shown in black contours in Figure 20 (left). It is possibly connected to 
the old pulsar PSR J1833-0827 (Gaensler & Johnston 1995). This pulsar is energetic enough 
to power HESS J1834-087, but its distance of 24 arc minutes renders an association unhkely; 
there is also no extended PWN detected so far. 



5.3.14. HESSJ1837-069 

This source (Figure 21) has an elongated shape with a semi-major axis of 7±1 arc 
minutes and a semi-minor axis of 3±1 arc minutes. With a flux of 13.4% of the Crab 
flux above 200 GeV, the energy spectrum can be well fitted by a power-law of photon 
index 2.27±0.06 (see Figure 21 (top right)). When using a large 9^ cut of (0.4°)^, the 
normalisation of the power-law fit increases by 30%, suggesting that the assumption of a 
Gaussian brightness profile does not hold for this source. HESS J1837-069 coincides with 
the southern part of the diffuse hard X-ray complex G25. 5+0.0. This source, which is 
12 arc minutes across, was detected in the ASCA Galactic Plane survey (Bamba et al. 

2003) . The nature of this bright X-ray source is unclear but it may be an X-ray synchrotron 
emission dominated SNR such as SN 1006 or a PWN. The brightest feature in the ASCA 
map (AX J1838. 0-0655), located south of G25. 5-1-0.0, coincides with the VHE emission. 
This still unidentified source was also serendipitously detected by BeppoSAX and also in the 
hard X-ray (20-100 keV) band in the Galactic Plane survey of INTEGRAL (Mahzia et al. 

2004) . 



6. Discussion 

Possible multi-wavelength associations of the fourteen new VHE 7-ray sources were 
discussed in the previous section and are summarised in Table 5. A further qualification 
from tentative associations to unambiguous counterpart identifications requires spatial or 
morphological coincidence, a viable 7-ray emission mechanism for the object, and consistent 
multi-wavelength behaviour according to the properties of the suggested identification. Most 
of the possible associations given above can be considered as plausible 7-ray emitters, but 
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only in a minority of these cases does a satisfactory positional coincidence exist. Four of the 
new sources apparently line up with SNRs. These SNRs have been detected and studied 
mostly in the radio domain, but in the cases of HESS J1640-465 and HESS J1813-178 X-ray 
emission was detected by ASCA, and there is marginal evidence for X-ray emission in ASCA 
archival data in the case of HESS J 1713-381. Only in the case of HESS J 1834-087 is there 
currently no evidence for associated X-ray emission. 

For three of the new sources, namely HESS J1616-509, HESS J1804-216 and HESS J1825- 
137, a sufficiently energetic nearby pulsar could account for the VHE 7-ray emission but 
the offset of the sources from the pulsar positions makes these associations less &m. Three 
of the new VHE 7-ray sources are consistent with the 95% positional error of unidenti- 
fied EGRET sources: HESS J 1640-465, HESS J 1745-303 and HESS J1825-137. Two of the 
new sources are located very close to a new class of highly absorbed hard X-ray sources 
detected recently by INTEGRAL (HESS J1632-478 and HESS J 1634-472). A third source 
(HESS J1837-067) may be related to X-ray emission of unknown origin detected by ASCA 
and INTEGRAL, which may belong to the same class of objects. The three remaining 
sources (HESS J1614-518, HESS J1702-420 and HESS J1708-410) have no plausible SNR, 
pulsar or EGRET counterpart and may belong to a new source class. Significant fiuxes of 
VHE 7-rays without accompanying X-ray and radio emission would suggest the absence 
of relativistic electrons and the presence of energetic nucleons. That such a population ex- 
ists is already being suggested by the detection of apparently "dark" accelerators by the 
HEGRA (Aharonian et al. 2002a) and H.E.S.S. (Aharonian et al. 2005b) collaborations. 
However, the possibility remains that these new sources are associated with SNRs or PWN as 
yet undiscovered at other wavelengths. Further multi- wavelength investigations are needed 
to clarify this picture. 

The nearest energy band relative to the H.E.S.S. energy range was covered by observa- 
tions with EGRET aboard the Compton Gamma-Ray Observatory. Its energy coverage was 
primarily between 30 MeV and 10 GeV and 7-ray source detections above 100 MeV resulted 
in the third EGRET (3EG) catalogue (Hartman et al. 1999). The only Galactic sources 
detected above 10 GeV are pulsars, for which the sparse photons at these energies correlate 
with the location and the pulse profile at lower EGRET energies (Thompson, Bertsch & 
O'Neal 2005). The bulk of photons above 10 GeV have to be attributed to Galactic (and 
extra-galactic) diffuse emission. 

When comparing the EGRET sources with the H.E.S.S. Galactic Plane survey result, 
it should be noted that even if an EGRET source spectrum is measured up to 10 GeV, more 
than a decade in energy is left uncovered below the H.E.S.S. energy threshold. Indeed, only 
a minor fraction of the H.E.S.S. sources coincide within the considerably larger location 
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uncertainty contours of EGRET sources. These sources are HESS J1640-465, HESS J1745- 

303, and HESS J1825-137, which arc consistent with the 95% location uncertainty of the 
unidentified EGRET sources 3EG J1639-4702, 3EG J1744-3011, and 3EG J1826-1302. In 
total 17 EGRET sources are catalogued within the H.E.S.S. survey region. Taking into 
account the 95% positional error contours, they cover an area of 9.76 deg^, resulting in a 
chance positional coincidence of a given H.E.S.S. source with an EGRET source of 32%. A 
correlated Galactic latitude distribution of the EGRET and H.E.S.S. sources, concentrated 
near the Galactic plane, could increase this probability. 

Physical arguments suggest that the EGRET source population and the set of new 
H.E.S.S. sources might be of different kinds. While the brightest 7-ray pulsars can be di- 
rectly traced up to ~ 20 GeV, EGRET data indicate the existence of GeV-cutoffs already in 
the 2-10 GeV range (Reimer & Bertsch 2001). The brightest pulsars, and several unidenti- 
fied EGRET sources with pulsar-like characteristics deviate from single power-law spectral 
representations (Bertsch et al. 2000), and are significantly better fitted when introducing 
high-energy spectral cutoffs. GeV-cutoffs arc indeed predicted from multi-frequency mod- 
elling of the 7-ray emission of pulsars (Cheng 2004; Baring 2004). Therefore, the emission of 
these unidentified EGRET sources may not extend up to the H.E.S.S. energy range. From 
Crab measurements by CELESTE (de Naurois et al. 2002), consecutively continuing the 
EGRET GeV-data, a characteristic change in the emission properties of the pulsar became 
apparent: pulsed emission is diminished or entirely absent, leaving only a steady emission 
component, in case of Crab attributed to the PWN. If the majority of EGRET sources in 
the Galactic Plane may be identifiable with 7-ray pulsars, this population will overlap with 
the H.E.S.S. energies only if a plerionic emission components exist. 

For PWN there are theoretical grounds to expect in some cases a displacement of the 
7-ray emission from the pulsar powering the nebula (Aharonian, Atoyan & Kifune 1997; 
Blondin, Chevalier & Frierson 2001). In such cases the association becomes more difficult 
to establish. Considering only pulsars with sufficient spin-down luminosity to explain the 
measured VHE 7-ray flux, four of the new sources have possible associations with pulsar 
wind nebulae. However, only in the case of HESS J1825-137, is there so far a convincing 
multi-wavelength argument for the offset from the pulsar position in the sense that an X-ray 
PWN offset in the source direction has been detected. 

The sensitivity shown in Figure 3 is that for point-like (< 0.1°) sources. The sensitivity 
of H.E.S.S. decreases with increasing source size: Fmin ~ -^point -source \/<^sourcc + '^psf/'^psf, 
where ctpsf ~ 0.08° for the analysis described here. As a consequence there is a bias in this 
survey against very extended sources. In particular, for sources larger than the size of the ring 
used for background estimation (0.6° radius) the flux threshold increases dramatically since 
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the off-source region used for the background estimation then also contains flux from the 
source. This effect also introduces a bias against close-by sources, potentially resulting in a 
latitude distribution narrower than that of the parent population. An additional consequence 
of the ring background technique is that this survey is insensitive to weak diffuse emission. 
A search for diffuse emission in this dataset will be reported elsewhere. Another effect that 
has an impact on the width of the latitude distribution is source confusion. Assuming that 
all sources he on the Galactic plane and using the mean RMS size of 0.1° of the new sources, 
the chance probability of an overlap between any two of the 17 sources in the scan region is 
estimated to be 6%. Therefore, whilst source confusion could in principle widen the observed 
latitude distribution, we expect this effect to be small. 

It is apparent from Table 5 that the most likely associations of the 7-ray sources lie at 
rather large distances (i.e. 4-10 kpc) within our Galaxy, and exhibit luminosities in the range 
3-20 xlO^^ erg s~^. Three of the new sources are located in a region of the Galactic plane 
that lines up with the 3 kpc arm tangent region (HESS J1640-465, HESS J1632-478 and 
HESS J1634-472), which hints at a distance of 8 kpc. HESS J1614-518 and HESS J1616-508 
are at a position in the plane that lines up with the Norma spiral arm. This arm appears to 
be one of the most important massive star forming regions of the Galaxy (Bronfman et al. 
2000). 

As it is clear from Table 5 it is difficult to explain the new VHE sources as belonging to a 
single source population. It is therefore not straight-forward to interpret the global properties 
of the sources. However, conventional Galactic particle accelerators such as young pulsars 
and SNRs have similar spatial distribution within the Galaxy. Molecular material and hence 
regions of star formation cluster rather close to the Galactic Plane, with an exponential 
scale height of <100 pc (Dame & Thaddeus 1994). 7-ray sources associated with stellar 
death (and birth) (SNRs, young pulsars, OB associations) should therefore follow similar 
distributions (see Figure 22). To compare such candidate populations with the observed 
Galactic latitude distribution of the survey sources assumptions must be made about the 
luminosity distribution, the radial distribution in the Galaxy, and the intrinsic size of the 
sources. A simple Monte-Carlo simulation shows that the measured b distribution of sources 
is rather insensitive to the assumed luminosity distribution. 

Taking the example case of SNRs as being a single class of counterparts to the new 
sources, it is possible to derive the expected properties of the detected 7-ray sources. For 
this exercise, the prescription of Drury, Aharonian & Volk (1994) as given in equation 1 is 
followed for the 7-ray emission of SNRs assuming supernova explosions with fixed energy 
lO^-*^ erg s"^ occurring at a rate of 1/40 year"-*^ and with a 7-ray emitting lifetime of 10^ 
years. In this case, 250 7-ray sources are expected within the Galactic disk. Assuming 
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the radial distribution of SNRs suggested by Case & Bhattacharya (1996), the measured 
number of the VHE 7-ray sources is consistent with expectations for plausible values of 
conversion efficiency into accelerated particles G = 0.2 and ISM density n = 1 cm~^. As the 
sensitivity of the survey is reduced for extended sources, some assumption of the intrinsic 
size distribution of sources must be made. Here, the Sedov solution is assumed to determine 
the size of the 7-ray emission region for SNR, of randomly sampled age, expanding into an 
ISM of density n. Reasonable agreement with the number and angular size of the detected 
sources are found for 1 < n < 10 cm~^. Figure 23 shows the expected distributions from this 
simple model for a scale height in the Galactic disk of 30 pc (red curve) and 300 pc (blue 
curve). The green curve is as the red curve but with 4 times as many sources, each with 
half the luminosity, illustrating the relative insensitivity of the shape of the distribution on 
the assumed source luminosity distribution. The difference between the blue and the red 
curve illustrates that the H.E.S.S. survey has a broad enough sensitive coverage in b to 
distinguish these two scenarios. For the 300 pc curve (blue), an intrinsic luminosity twice 
that of the 30 pc case has been assumed to give reasonable agreement with the number of 
detected sources. As can be seen from Figure 23 the experimental data favour a scale height 
for the parent population of < 100 pc in line with expectations for young SNRs and PWN 
or massive star forming regions. 



7. Conclusions 

The inner part of our Galaxy now contains seventeen known sources of VHE7-rays, 
compared to three before the H.E.S.S. survey. The increased number of sources allows us 
to consider the behaviour of populations of sources for the first time in this wave-band. A 
major task of multi-wavelength follow-up observations and archival searches now lies ahead to 
understand the processes at work in these astrophysical particle accelerators. The paradigm 
of cosmic-ray acceleration in SNRs is consistent with our new findings but it seems clear 
that the new sources arc not drawn from a single population. Follow-up observations of the 
new sources with H.E.S.S. are planned as well as an extension of this survey into unexplored 
regions of the Galactic Plane. 
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Table 5: Summary of possible counterparts to the H.E.S.S. VHE 7-ray sources in the inner 
Galaxy. The ofTset denotes the angular distance of the VHE 7-ray source from the possible 
counterpart. Distances are taken from Green (2004); Manchester et al. (2005) and references 
therein. The last column gives the imphed luminosity between 0.2 and 10 TeV. The three 
sources that were known before the H.E.S.S. survey in the inner Galaxy (HESS J1747-218, 
HESS J1745-290 and HESS J1713-397 have been added for completeness to this table. The 
possible source classes have been abbreviated according to the following scheme: Supernova 
remnant (SNR), unidentified sources (UID), pulsar wind nebula (PWN), black hole candidate 
(BH), X-ray binaries (XRB). 



Name 


Possible 


Class 


Offset 


Distance 


Luminosity 




Counterpart 




(arcmin) 


(kpc) 


(10^^ erg s"^) 


J1614-518 












J1616-508 


PSRJ1617-5055 


PWN 


10.4 


6.5 


20.2 


J 1632-478 


1GRJ16320-4751 


XRB 


3 






J 1634-472 


IGRJ16358-4726 / G337.2+0.1 


XRB/SNR 


13/10 


-/14 


-/28.3 


J1640-465 


G 338.3-0.0/3EG J1639-4702 


SNR/UID 


0/34 


8.6 


16.2 


J1702-420 












J1708-410 












J1713-381 


G 348.7+0.3 


SNR 





10.2 


5.2 


J1713-397 


RXJl 713. 7-3946 (G 347.3-0.5) 


SNR 





1 


1.2 


J 1745-290 


Sgr A East/Sgr A* 


SNR/BH 





8.5 


12.8 


J 1745-303 


3EG J 1744-30 11 


UID 


10 






J1747-281 


GO. 9+0.1 


PWN 





8.5 


4.4 


J1804-216 


G8.7-0.1/PSR J1803-2137 


SNR/PWN 


21/10.8 


6/3.9 


16.5/7.0 


J1813-178 


G 12.82-0.02 


SNR 





4 


3.1 


J1825-137 


PSR J1826-1334/3EG J1826-1302 


PWN/UID 


11/43 


3.9/- 


6.1/- 


J1834-087 


G 23.3-0.3 


SNR 





4.8 


4.4 


J1837-069 


AXJ1838.0-0655 


UID 


6 
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Fig. 14. — Smoothed excess map (left, smoothing radius 0.05°) of the region surrounding 
HESS 1708-410 along with a nearby SNR which was considered as counterpart. The differ- 
ential energy spectrum (top) and the radial distribution (^^) plot (bottom) are shown on the 
right hand side. 
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Fig. 15. — Smoothed excess map (left, smoothing radius 0.07°) of the region surrounding 
HESS J1713-381 along with nearby SNRs which were considered as counterparts as discussed 
in the text. In the bottom right corner of this excess map, the strong VHE emission from 
SNR RX J1713. 7-3946 is apparent. The white contours show the radio emission as seen 
in the Molonglo Plane survey at 843 MHz (Green et al. 1999). The marginal excess in 
positional agreement with G 348.5+0.1 (CTB 37 A) is not significant after accounting for 
all trials for the analysis presented here. The differential energy spectrum (top) and the 
radial distribution (6'^) plot (bottom) are shown on the right hand side. Also in the radial 
distribution, some indication for an excess ~ 0.4° away from HESS J1713-381 can be seen, 
due to the marginal excess at the position of CTB 37 A. The dashed white contours show 
the 95% positional error of the unidentified EGRET source 3EG J1714-3857. 



- 49 - 




B= (deg^) 



Fig. 16. — Smoothed excess map (left, smoothing radius 0.10°) of the region surrounding 
HESS J1745-303 along with nearby pulsars and SNRs which were considered as counterparts 
as discussed in the text. The differential energy spectrum (top) and the radial distribution 
{6^) plot (bottom) are shown on the right hand side. 




Fig. 17. — Smoothed excess map (left, smoothing radius 0.06°) of the region surrounding 
HESS J1804-216 along with nearby pulsars and SNRs which were considered as counterparts 
as discussed in the text. The white contours show the X-ray emission from G8. 7-0.1 as 
detected by the ROSAT satellite, the black contours show the 20 cm radio emission detected 
by the VLA. The differential energy spectrum (top) and the radial distribution (6*^) plot 
(bottom) are shown on the right hand side. 
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Fig. 18. — Smoothed excess map (left, smoothing radius 0.04°) of the region surrounding 
HESS J1813-178 along with nearby pulsars and SNRs which were considered as counterparts 
as discussed in the text. The white contours shows the smoothed X-ray count map for the 
region surrounding AXJ1813-178 obtained by ASCA satellite with the GIS instrument in 
the 0.7-10 keV band, the black contours show 20 cm radio emission as observed by the VLA. 
The strong radio source to negative galactic longitudes from HESS Jl 8 13-178 is connected to 
the W 33 complex and extends further than shown here. The excess seen to negative Galactic 
longitudes at the edge of the ASCA GIS field of view is likely due to contamination from 
GX13+1. The radial (lower right) distribution {6"^) plot shows that the source is marginally 
extended with respect to the PSF of the H.E.S.S. system. 
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Fig. 19. — Smoothed excess map (left, smoothing radius 0.09°) of the region surrounding 
HESS J1825-137 along with the nearby pulsar which was considered as counterpart as dis- 
cussed in the text. The differential energy spectrum (top) and the radial distribution (6*^) 
plot (bottom) are shown on the right hand side. 
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Fig. 20. — Smoothed excess map (left, smoothing radius 0.08°) of the region surrounding 
HESS J1834-087 along with nearby pulsars and SNRs which were considered as counterparts 
as discussed in the text. The positional coincidence with SNR G23.3-0.3 (W41) is evident, 
the black contours show 20 cm radio emission recorded by the VLA (White, Becker & Helfand 
2005). The differential energy spectrum (top) and the radial distribution {O"^) plot (bottom) 
are shown on the right hand side. 
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Fig. 21. — Smoothed excess map (left, smoothing radius 0.05°) of the region surrounding 
HESS J1837-069 along with nearby SNRs which were considered as counterparts as discussed 
in the text. The black contours shows the smoothed X-ray count map for the region sur- 
rounding AX J1838. 0-0655 obtained by the ASCA satellite with the GIS instrument in the 
0.7-10 keV band. The positional coincidence of HESS J1836-069 with the ASCA source 
AX J1838. 0-0655 is evident. The differential energy spectrum (top) and the radial distribu- 
tion (0^) plot (bottom) are shown on the right hand side. 
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Fig. 22. — Distribution of Galactic latitude of the seventeen VHE 7-ray sources as given in 
Table 5 (the fourteen new plus the three previously known). The mean of the distribution is 
-0.17 with an RMS of 0.31. On the left, the distribution is shown along with the distribution 
in Galactic latitude of possible counterpart populations to the VHE 7-ray sources. SNRs 
from Green (2004) are shown in dash-dotted purple, high spin-down luminosity pulsars 
{E > 10^^ erg s~^) taken from the ATNF catalogue (Manchester et al. 2005) are shown as a 
dashed blue curve. On the right the distribution of VHE 7-ray sources as on the left is shown 
along with the latitude distributions of matter in our Galaxy. The distribution of molecular 
hydrogen (HII) (Paladini et al. 2003) is shown in dash-dotted red, the distribution of atomic 
hydrogen (HI) (Dickey & Lockman 1990) is shown in dotted purple and the distribution of 
CO (Dame, Hartman & Thaddeus 2001) is shown as dashed green line. 
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Galactic Latitude (deg.) 

Fig. 23. — Comparison of the Galactic latitude distribution of the VHE 7-ray sources in 
the survey region with the prediction of a simple source population model assuming scale 
heights in the Galactic disc of 30 pc (red and green curves) and 300 pc (blue curve). The 
sources are assumed to be mono-luminous. The curves are not normalised and give the 
predicted number of detected sources in each latitude bin taking into account the reduced 
sensitivity at higher latitudes 3. The legend gives the physical parameters in each case with 
Lq = 1.3 X 10^^ erg s~^ (in the energy range 0.2 TeV to 10 TeV) as previously defined. 
denotes the number of expected 7-ray sources. 



